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Cascades of Particles Moving at Finite Velocity 
in Hyperbolic Spaces 
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C^) Abstract 

A branching process of particles moving at finite velocity over the geodesic lines of the hyperbolic 
space (Poincare half-plane and Poincare disk) is examined. Each particle can split into two particles 

i— j only once at Poisson paced times and deviates orthogonally when splitted. At time t, after N(t) 

Poisson events, there are N(t) + 1 particles moving along different geodesic lines. We are able to 
obtain the exact expression of the mean hyperbolic distance of the center of mass of the cloud of 
particles. We derive such mean hyperbolic distance from two different and independent ways and we 
study the behavior of the relevant expression as t increases and for different values of the parameters 
c (hyperbolic velocity of motion) and A (rate of reproduction). The mean hyperbolic distance of each 

CLj moving particle is also examined and a useful representation, as the distance of a randomly stopped 

_J particle moving over the main geodesic line, is presented. 

Keywords: Branching processes, difference-differential equations, hyperbolic Brownian motion, hyperbolic 
Cj trigonometry, Laplace transforms, non-Euclidean geometry, random motions. 
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O^ 1 Introduction 



Random motions in hyperbolic spaces have been studied since the Fifties and much emphasis has been 
placed on the so-called hyperbolic Brownian motion on the Poincare half-plane (see, e.g., Gertsenshtein 
and Vasiliev [3J, Getoor [3], Gruet |B], and Lao and Orsingher [5]). 

Hyperbolic Brownian motion has been revitalized by mathematical finance since some exotic financial 
products (Asian options) have a strict connection with the stochastic representation of the hyperbolic 
Brownian motion (Yor [Ll]). 

Branching hyperbolic Brownian motion has been analyzed by Lalley and Sellke [10] who investigated 
the connection between the birth rate and the underlying dynamics in supercritical and subcritical cases. 
Also Kelbert and Suhov [7J , [5] have studied the asymptotic behavior of the hyperbolic branching Brow- 
nian motion, developing the ideas in [3J and [TO] . 

The space on which the above considered hyperbolic Brownian motions develop is the Poincare half- 
plane (and its higher-dimensional equivalents, see Gruet [5]) or the Klein model (see [7J and [B]). 

The half-plane Poincare model is a fine tool to describe the light propagation in a non-homogeneous 
medium where, on the basis of Fermat's principle, the angle a(y) between the tangent to the geodesic curve 
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at point of ordinate y satisfies the equality [sin a (?/)]/[«/] = 1/k. The randomly scattered irregularities 
in the medium cause deviations of the trajectories and their decomposition into different rays. 

Random motions with finite velocity have been considered in Orsingher and De Gregorio |12j (on H^ 
and on the Poincare disk) with the assumption that they develop on geodetic lines and have independent 
components. 

Several models of random motions in H^ with finite velocity have been examined in Cammarota and 
Orsingher [lj, where the components of the motion have been assumed dependent and the particle moves 
on mutually orthogonal geodesic lines. 

Here we study a random motion of a cloud of particles moving at finite velocity on geodesic lines of 
the hyperbolic space H7£ . Such particles are generated by successive splitting out of a unit-mass particle 
initially placed at the origin O of H% . The disintegration process of the unit-mass particle is governed 
by an underlying Poisson process of rate A as follows. The original particle keeps moving on the main 
geodesic line with constant hyperbolic velocity c; at the first Poisson event it breaks into two parts of 
equal mass 1/2. One particle continues its motion on the main geodesic line, while the other one deviates 
orthogonally. At the second Poisson event the deviated particle is separated into two pieces each of mass 
1/2 2 ; one piece continues its motion on the same geodesic line whereas the other one starts moving on 
the geodesic line orthogonal to that joining its position with the origin O. In general (see Figure [l]), at 
the fe-th Poisson event, only the deviating particle of mass l/2 fc breaks into two fragments of equal mass 
l/2 fc+1 ; the first continues its motion on the same geodesic line while the other one deviates orthogonally. 

If up to time t, N(t) Poisson events have occurred, we have N(i) + 1 particles running, at a constant 
hyperbolic velocity c, along different geodesic lines with a mass depending on the instant of separation 
from the generating particle. 

In the above branching process each particle can reproduce only once and the particle splitting at 
time Tfc of the fc-th Poisson event and which never more disintegrates will preserve its mass, equal to 
l/2 fc , for the successive time interval (Tk,t). 

Our main result concerns the dynamics of the center of mass cm of the cloud of particles performing 
the branching and diffusion process. In particular, we are able to give an exact expression for the mean 
hyperbolic distance from the origin O, r] cm (t), of the center of mass cm at any time t > 

Mcoshwft)} = { , + , } (1.1) 
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We give two different and independent proofs of the above result: our first technique is based on Laplace 
transforms, while the other one brings about the following non-homogeneous second-order differential 
equation 
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which is satisfied by (1.1) 



The behavior of the hyperbolic distance of each individual particle can be compared with result ( 1.1 ). 
In a previous paper (see pQ) we have shown that the mean hyperbolic distance r)(t) of the particle which 
underwent changes of direction at all Poisson events is 
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E{cosh V (t)} = ; Ce J-g J + , 6 . (1.2) 
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If in ( 1.2 1 the Poisson rate A is replaced by A/2 we see that the exponential terms in ( 1.1 ) and ( 1.2 ) have 
the same form but different weights. 

We also examine the mean hyperbolic distance of each individual particle which stops changing direc- 
tion after the fc-th Poisson event. Our main result shows that 

1 



E {cosh r] k (t) I {N{t )> k} } = -£ / coshc(£ - s)h(k,c,s)g(s; k, A) ds (1.3) 
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where g(s;k, A) = e "]d~J~ is a Gamma distribution and h(k,c,s) = J2 r =o (r)^ Y r k {e cs ( 2Yr - k ^} with 
Y Tt k ~ Beta(r, k — r) so that (2F r ,fc — 1) € (— 1, 1) and with the assumption that Yq^ — 1 and 5^fc,fe = — 1. 



Thus ( 1.3 1 shows that the mean hyperbolic distance, at time t, of a particle generated at the fc-th Poisson 



event can be seen as the mean hyperbolic distance of a particle which never deviates from the main 
geodesic line and which starts moving at a random time with law h(k, c, s)g(s; k, A). 



2 Some geometrical features of the hyperbolic spaces 

We present in this section some basic features of the Poincare half-plane H^ = {(x, y) : y > 0} which is 
endowed with the metric 

ds= V(d^ + W 
y 

Some informations on hyperbolic spaces and non-Euclidean geometry can be found in Faber [5] and 
Meschkowski [llj . The position of points in H^ can be given in Cartesian coordinates (x,y) or in 
hyperbolic coordinates (j], a). These are connected by means of the well-known relationships 



x 



sinh rj cos a 



cosh rj — sinh rj sin a ' 
" cosh r\— sinh rj sin a ' 



7?>0, 



f <a<f. 
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The hyperbolic coordinate 77 represents the distance of (x, y) from the origin O = (0, 1) of H^ measured 
by means of the metric (2.1 1. The coordinate a is the angle of the tangent in O to the half- circumference 



joining O with (x, y) (see Rogers and Williams [T3] page 213 and Cammarota and Orsinghcr [1 for some 

details). 

For each point (x, y) £ H^ it is possible to obtain the hyperbolic distance 77 as well as the angle a by 

means of the formulas 



cosh 77 = 



x 2 + y 2 
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Formulas (2.3) are easily derived from (2.2) (see [H]). We can obtain formulas (2.2) from (2.3) as follows. 
By substituting 



x = tana ± v/tan a + 1 — y 2 , 



in the first relationship of (2.3) we have that 



y cosh 77 — 1 — tan a = ± tanaytan 2 a + 1 — y 2 , 
and after some manipulations we arrive at 
= (y cosh 77 — 1)" — y sin a sinh 77 = (y cosh 77 — 1 — y sin a sinh 77) (77 cosh 77 — 1 + y sin a sinh 77), 



which yields the second formula of (2.2). Since from (2.3) y cosh 77 = a; tan a— 1, also the second relation- 



ship of (2.2) immediately follows. 



In HJ the trigonometrical formulas we need are either the Pythagorean theorem for right triangles 



cosh 77 = cosh 77! cosh 772, 
or its Carnot extension for arbitrary triangles 

cosh 77 = cosh 771 cosh 772 — sinh 771 sinh 772 cos(ai — 02). 



(2.4) 
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Clearly, if a.\ — a.2 = tt/2, formula (2.5) becomes (2.4) 



The half-plane H% can be mapped onto the disk D = {(u, v) 



mapping 
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-v 2 < 1} by means of the conformal 
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The rr-axis of H^ i s mapped onto the boundary dD of D while the origin O is transformed into the center 
of D. An arbitrary point (x, y) £ H^ is mapped onto the point (u, v) £ D with coordinates 

2x x 2 + y 2 — 1 



x 2 + (y + l) 2 ' x 2 + (y + l) 2 



A point (u, v) £ D is instead mapped by (2.6) into the point (x, y) with coordinates 



_ 2w _ 1 - (u 2 +v 2 ) 

X ~ u 2 + (1 - v) 2 ' V ~ u 2 + {l-v) 2 ' 

A similar mapping is the so-called Cayley transformation which reads 



w= t — Z , (2.7) 

i + z 



and it slightly differs from (2.6 ). The geodesic lines of H^ with radius r and center at (xq, 0), are mapped 



by (2.6) into arcs of circumferences inside D with center at 

2xq Xn — r 2 — 1 



and with radius R given by R 2 = 



•o 

4r 2 
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3 Description of the randomly moving and branching model 

We assume that a unit-mass particle is placed at time t = at the origin O of H2 and starts moving on 
the main geodesic line represented in H^ by the half-circle of radius 1 passing through O. This particle 
chooses with probability 1/2 one of the two possible directions and moves with constant hyperbolic 
velocity equal to c (see Figure HI (a)). The hyperbolic velocity 




is assumed to be constant. For an Euclidean observer, the closer to the a;-axis is the moving particle the 
slower it moves. 

A Poisson process of rate A governs the changes of direction. At the first Poisson event the particle 
splits into two pieces of equal mass: one continues its motion on the same geodesic line while the other one 
starts moving (in one of the two possible directions) on the geodesic line orthogonal to the previous one 
(see Figure [ll (b)). In general, at the fc-th Poisson event, the deviating particle of mass l/2 fc undergoes 
a further decomposition: one splinter of mass l/2 fc+1 continues undisturbed its motion, while the other 
one, also of mass l/2 fc+1 , is forced to move onto the geodetic line orthogonal to that joining O with the 
position it occupied at the time where the splitting took place. 

Therefore, if no Poisson event occurs (i.e., {N(t) — 0}, where N(t) is the number of Poisson events 
in [0, £]), the unit-mass particle is located, at time t, on the first geodesic line at an hyperbolic distance 
from O equal to T]o(t) — ct (see Figure fl] (a)). 

If one Poisson event happens, at time S\ < t (i.e., {N(t) — 1}), then, at time t, one fragment of mass 
1/2 will be at distance rjo(t) — ct on the first geodesic line, while the other splinter, also of mass 1/2, will 
be located at hyperbolic distance r)i(t) given by 

cosh 771 (£) = cosh c Si cosh c(t — Si ) , (3.1) 



where in formula (3.1) the hyperbolic Pythagorean theorem (2.4) has been applied (see Figure fl] (b)) 



If N(t) = n, the process described above produces n + 1 splinters. The particle generated at the fc-th 
Poisson event, k = 0, • • • n — 1, (and which will not break up after the fc-th event) has mass l/2 fc+1 and 
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Figure 1: In (a), the trajectory of the unit-mass particle initially placed at the origin O of -H^" and moving on 
the main geodesic line is shown. When N(t) = 0, no disintegration occurs. In (b), (c), (d), and (e), the 
trajectories of the particles generated by N(t) = 1,2,3,4 Poisson events are plotted. The relevant mass 

associated with each particle is also indicated by a suitable label. In (f), for each particle the relevant hyperbolic 
cosine of the hyperbolic distance from the origin is also reported. 



is located, at time t, at the hyperbolic distance r)k(t) from O. Such a distance, in force of the hyperbolic 

Pythagorean theorem, reads 

fc+i 

coshryfc(t) = TT coshc(Sj — Sj-i), 
i=i 

where So = 0, Sk+i = t and the random times Sj with j = 1, • • • k represent the instants where the 
Poisson events happen and the deviations of motion occur. The last splinter, which has changed direction 
at all fission events, has mass 1/2™ and, at time t, is at an hyperbolic distance from O given by 

n+l 



coshr] n (t) = TT coshc(Sj — Sj-i). 



where So = and S„ + i = t (see Figure [I] (f) ) . 

In general, if the number of splits recorded is N(i), the number of particles is N(t) + 1 and each runs 
on a different geodesic line of H^- The hyperbolic distance of the center of mass cm of the cloud of 
moving particles, at time t > 0, is denoted by r] cm (t) and is represented by 

N(t)-1 _, fe+l _, N(t) + 1 



coshr] cm (t) = Y^ 2k+i II coshc ( S j ^ S j-i) I {N(t)>o} + ^^y J| coshc ( S i - S i-i)> 



(3.2) 
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where Si, S2 . . . Sjv(t) are the random times at which the disintegrations occur, and So = 0, Sjv(t)+i = t. 



The second term in (3.2 1 refers to the splinter which underwent all disintegrations occurred until time t, 
while the first one is related to those particles produced during the branching process. 

The assumption that the particles deviate on geodesic lines orthogonal to those joining the origin 
O with their current position is crucial since it makes the hyperbolic Pythagorean theorem applicable. 
Otherwise it should be necessary to apply the Carnot hyperbolic formula and this would make the analytic 
treatment of the problem extremely difficult. 

Under the condition that N(t) — n, the mean hyperbolic distance of the center of mass of the cloud 
of particles at time t is 

7i— 1 -. fe+1 -. n+1 

£{cosh?7 cm (t)|iV(i) = n] = J2 tf+i II cos M^ - S,-_i)J { „ >0} + — J] coshc(5, - S,_i). (3.3) 
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We observe that the n instants S\ , • • • ,S n where the fissions take place are uniformly distributed under 
the condition that N(i) = n, and possess density 

n! 
Pr{Sx € dsi, • • • ,S„eds n }= — dsi • • • ds n 

t n 

for = So < 8\ < ■ ■ ■ < s n+ i = t. Therefore, under the condition that the number of splitting events is 
N(t) = n, the hyperbolic distance r)k(t) of the fc-th splinter is for k = 0, • • • n — 1 

I pt pt pt K+l 

E{coshr) k (t)\N(t) = n} = -^ / dsi • • • / ds k • • • / ds n JJ cosh c(Sj - Sj_i) 

1 Jo Js k - ± Js n -X j =1 

= Wo dSl "J Sk J k (n-k)! n coshc ( s i- 5 ^) 

- ^G„, fe (t), (3.4) 



and for k = n 

E{coshr)„(t)\N(t) = n} 



n! /"* /"* " +1 

— / dsf- ds„ TT cosh c(s_j - Sj_i) 

= ^G„,„(t). (3.5) 

The branching process described above can be adapted to the Poincare disk _D: since the mappings 



(2.6 1 and (2.7) preserve the hyperbolic distance, the trajectories of splitting and moving particles can be 
conveniently depicted in D as in H^ (see Figure 111 and 0) . 

We restrict ourselves to the mean hyperbolic distance because this leads to fine explicit results. The 
analysis of the distribution of the hyperbolic distance, even in the case of a random motion of an individual 
particle that changes direction at all Poisson events, implies a much more complicated analysis and is 
almost intractable since multiple integrals of the form 

E{e iacoshT >M\N(t) = n} = — f ds x • • ■ f e ^n^h°^^ 3 -^-i) ( i Sn 

t U JO ij„_! 

must be evaluated (see Cammarota and Orsingher [1] for details on this point). 

4 Mean hyperbolic distance of the system of randomly moving 
and disintegrating particles: the Laplace transform approach 



We are able to obtain the explicit form of the mean-value of ( |3.2[ ) by means of two different and inde- 
pendent approaches. In this section we present the Laplace-transform derivation which leads to our first 
theorem. 
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Figure 2: Same trajectories as in Figure fil represented through the Poincare disk model. 



Theorem 4.1. The mean-value of the hyperbolic distance (3.2) of the center of mass is 



E{cosh.r] cm {t)} = 
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t > 0. 



Proof 



In view of (3.3 1, (3.4) and (3.5), wc first note that 



E {cosh Vcm (t)\N {t) = n} = 
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Our task is therefore to study the Laplace transform 
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where /i > 0. We evaluate the Laplace transform appearing in (4.2) in the following way. If 7 = A + \x 



and c < 7, then we have, by successively inverting the inner integrals, that 
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In the last step above the change of variable Sj — Sj-i — Wj applied k times leads to the final expression. 
For k — n, from the previous calculations, we obtain 
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we have that 
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(4.3) 



where the last sum converges if /j, satisfies the inequality 2c 2 < A 2 + 2fi 2 + 3A/i. The double sum in (4.3) 
can be calculated by inverting the order of summation in the following way: 
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The inversion of the Laplace transform is made possible by suitably rearranging the expression (4.4) as 
follows: 
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By means of the decomposition 
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the expression (4.5) can be further worked out by writing 
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2VA 2 + 2 4 c : 




A 





,, , 3A-VA 2 +2 4 e 2 ,, , 3AWA 2 + 2'c 2 

A* T 2 2 ' ' 2 2 



A + 2c A - 2c 



fi — c /i + c 



.. , 3A-VA 2 +2 4 c 2 , SA+TA^T 

A* T 2 2 ft 2 2 



2'r 2 



2 2 (2 M +|A)' 



(4.6) 



It is now a simple matter to invert the Laplace transform (4.6 1 and we arrive at the mean hyperbolic 
distance of the center of mass in an integral form 



E {cosh r) cm (t)} 

'1 A 

.2 + VA 2 +2 4 c 2 
A(A + 2c) 

h 2~7I 2 + 2 4 c 2 
A(A - 2c) 

^27T 2 + 2 4 c 2 



, 3X— v/X 2 +2 4 

„ — t Z7i 



- (t - s) M=Vpg 



e - cs e" ( ^ s) i 2 ds 



2 VA 2 + 2 4 c 2 

r* 

ds — I e~( 

Jo 

rt 



, 3A + Va 2 +2 4 



, ^ 3A + Va 2 +2 4 c 2 



,. n 3A + VA 2 +2 4 c 2 

e - cs e- {t - s) *" ds 



23 



^A,-c-») 3> -y 



ds 



_ 3 As _ (t _ s) 3A + VA 2 + 2lc2 
g 2 2 g V / 2 2 ,^ g 



(4.7) 



The expression (4.7) can be further developed and simplified by observing that, after some simple calcu- 



lations, we have that 

A(A + 2c) r '* 



-(t-s) ^-^ 



2y/X 2 + 2 4 c 2 [Jo 
A(A + 2c) 



2VA 2 + 2 4 c 2 2 2 c + 3A - VA 2 + 2 4 c 2 

A(A + 2c) 2? 

~ 2VA 2 + 2 4 c f 2 2 c + 3A + VA 2 + 2 4 c 2 



ds — i < ' < 



e ct - e- 1 



y ' 22 ds 



3A-y / A 2 +2 4 c 



e ct -e 



, 3A + Va 2 +2 4 c 2 
_t 72 



(4.8) 
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Similar manipulations yield 



A(A - 2c) 
2VA 2 + 2 4 c 2 



e-Ce-tM ^I^ 



ds — I < < 

/o 



CS„ It Sj 2 2 ^ 



A(A - 2c) 






2 2 




2VA 2 + 2 4 c 2 

A(A - 2c) 


-2 2 c 


+ 3A 


-VA 2 
2 2 


+ 2 4 c 2 [ 






3A-Va 2 + 2 4 c 2 



2VA 2 + 2 4 c 2 -2 2 c + 3A + VA 2 + 2 4 c 2 



e- ct - e~* 



3A + VA 2 | 2 ',- 



(4.9) 



and also 



A 



_3 As -((-s)^p3j , / -^As -(t-si- = --. | 

s 55 e y I 22 dS + / e 22 e *■ I 2= ds 



e 22 e 



A 
2VA 2 + 2 4 c 2 



. 3A + Va 2 +2 4 c 







-4rA< 

e 2 2 



, VA2 +2 4e2 ./g , ,, ! , 



(4.10) 



By inserting results (4.101, (4.8), and (4.9) into (4.7) we now obtain the final formula 



E {cosh r] cm (t)} 

-4rAt 

e 22 



3A 



VF+2V 



-WA2+2 4 c 2 



1- 



3A 



--WA 2 +2 4 e 2 

e 2^ v ^ 



A + 2c 



2A 



. 3A-VA 2 . _■ 1, 



2(A + 3c) e VA 2 + 2 4 c 26 



A -2c 



e- ct + 



2A 



, 3A + Va 2 +2 4 c 2 



2(A - 3c) VA 2 + 2 4 c 2 



VA 2 + 2 4 c 2 , 

A - 2c A + 2c 

2 2 c - 3A + \/A 2 + 2 4 c 2 2 2 c + 3A - VA 2 + 2 4 c 2 

A + 2c A - 2c 

2 2 c + 3A + VA 2 + 2 4 c 2 2 2 c - 3A - \A 2 + 2 4 c 2 



-4At -WA 2 +2 4 c 2 

e 22 e2 2 v T 



3A 



3A 2 - AVA 2 + 2 4 c 2 - 2 3 c 2 



2 2\/A 2 + 2 4 c 2 VA 2 + 2 4 c 2 V 5A - 3\/A 2 + 2 4 c 2 



-iA( --WA 2 +2 4 c 2 

+e 22 e 22 v T 



1 



3A 



3A 2 + A\/A 2 + 2 4 c 2 - 2 3 < 



2 2VA 2 + 2 4 c 2 VA 2 + 2 4 c 2 V 5A + 3 VA 2 + 2 4 c 2 



A + 2c ct A -2c _ ct 

e + -T-, -^e ct 



2(A + 3c) 2(A-3c) 

VA 2 + 2 4 c 2 



5A + 3\/A 2 + 2 4 c 2 5A - 3\/A 2 + 2 4 c 2 



A + 2c ct A -2c _ ct 

-e + — — e . 



2(A + 3c) 2(A-3c) 



Remark 4.1 . Apparently a critical point in formula (4.1 1 is A = 3c. We show that the mean hyperbolic 



distance (4.1 ) is finite for A = 3c. We first write e = A — 3c and r(e) = \Je 2 + 6ce + 5 2 c 2 and evaluate 
the limit 



A -2c 



2 3 c 2 e ^ ; 



A™c2(A-3cf VA 2 + 2 4 c 2 5A - 3%/A 2 + 2 4 c 2 

e + c _ ct 2 3 c 2 

-e 



lim 

e^o 2e 



r(15c + 5e — 3r) 



e -^(9c+3£-r) 



(4.11) 



which refers to the components of (4.1 ) with diverging coefficients. We expand the second exponential in 



( 4.11[ ) as 

e -5r(9c+3e-r) 



e -ci e -4r(5c+3e-r) _ & -ct 



t t 2 

1 - 7^(5c + 3e - r) + -e(5c + 3e - r) 2 + o(e 3 ) 
2 Z 2° 
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Figure 3: In (a) and (b) the mean-value _E{cosh?7 cm (i)} is plotted for c = 1 and different values of A, the curve 

corresponding to A = 3c is obtained by plotting the limit value 

£{cosh7y cm (t)} = e" ct (^ffj + ^ct) + e^ + e" 2 ? J^. 



where we have taken into account that r(e) 
form 



5c for e — > 0. It is now convenient to write (4.11 ) in the 



lim 



e + c 

2e 



r(15c 

By considering that r{e) ~ 5c 



2 3 c 2 

+ 5e- 



3r) 



2£c 2 (5c + 3e-r) 
r(15c + 5e — 3?') 



i 2 c 2 (5c + 3e 



lim 

e^0 



2e r(15c + 5e - 3r) 



|e for e 

_ 2-7 
"15 2 "' 



2 2 r(15c + 5e-3r) 
0, we can easily realize that 



+ o(s 2 ) 



2fc 2 (5c + 3e-r) _ 3 
e-s-o r(15c + 5e - 3r) ~ 2 ■ 5 



ci, lim 



t 2 c 2 (5c + 3e~r) 2 
o 2 2 r(15c + 5e - 3r) 



This suffices to show that near A = 3c the mean hyperbolic distance is finite. Considering also the two 



additional terms of (4.1 1 leads to the following asymptotic estimate of the mean hyperbolic distance of 



the center of mass for large values of t and near A ~ 3c. In particular, we have that 

5 
E{coshr) cm {t)} ~ — e c *. 

In Figure|3j the function E {cosh r] cm (t)} is plotted for c = 1 and different values of A in the neighborhood 
of A = 3c (including the limiting case A = 3c). 

Remark 4.2. We now show that the center of mass (as well as each individual particle) goes further and 
further away from the origin O of H^ (or equivalently it migrates towards the frontier of the Poincarc 
disc D) as time t increases. In order to show this result we work on (4.7) and, after some manipulations, 
we obtain that 



d* 



E{cosh.r} cm {t)} = 



2 2 cV 



-A( 



t 



sinh — vA 2 

2 4 c 2 2 2 V 



2 4 c 2 



2cA 2 



VF+2¥ 

2 2 c 2 A 
VA 2 + 2 4 c 2 



e ^ sinh c(t — s) sinh —^ v A 2 



2 4 c 2 di 



/ e~^ Xs cosh c(t - s) sinh -i \/\ 2 + 2 4 c 2 ds. (4.12) 

Jo 2 



The first term in the right hand side of (4.12| is produced by the derivatives of all terms of (4.7) 



while the integrals stem from the third and fourth term only. From (4.12) it is easy to check that 



^2 E {cosh r] cm (t)} > for all t so that the center of mass gets off from the starting point with positive 
acceleration. 
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c=0.05, t=40, N(t}=20 



c=0.05, t=40, N(t)=20 



c=0.05, t=40. N(t)=20 






3.05, t=50, N(t)=300 



3.05, t=50, N(t)=300 



3.05, t=50, N(t) = 300 



Figure 4: In (a), (b), and (c) the trajectories and the positions of the splinters at time t = 40, when c = 0.05 

and N(t) = 20, are shown. In (d), (e) and (f) only the positions of the splinters at time t = 50, when c = 0.05 

and N(t) = 300, are drawn. In (a) and (d) each splinter chooses the clockwise direction, in (b) and (e) each 

splinter chooses the counterclockwise direction, in (c) and (f) the clockwise and counterclockwise directions are 

alternatively chosen. 

5 Equation governing the hyperbolic distance 



We are able to confirm result (4.1 1 by a completely different method based on a system of differential 



equations governing all the quantities appearing in the mean hyperbolic distance of the center of mass 



oo n—1 

E{cosh Vcm {t)} = e- xt V A" ^ 



1 



71=1 



fc=0 



G n ,k{t) + ' 



-\t v^ 



A" 



^ 2' 

ra=0 



(*n,n{t)- 



(5.1) 



We start by deriving the difference-differential equations governing the functions G n ^(t), with t > and 
< k < n. As it has been pointed out in our previous paper [T], the j'-th moment of the conditional 
hyperbolic distance of a single particle leads to a difference-differential equation of order equal to j + 1. 
The same type of phenomenon occurs in the present case where the hyperbolic distance of a cloud of 
points is envisaged. 



Lemma 5.1. The functions represented by the following multiple integrals 



G n ,k(t) 



d.s 



l • • ■ 



dsi 



(t-Bk)' 



-k fe 

FT cosh ( 



(n-k)\ 
are solutions to the difference- differential equations 



s i-i 



) coshc(i — Sfc), 



< k < n, 



di^Gn.k — 2gjG„_i i fc — G n -2,k + C Lrn,fe) 
d 2 

d/- 



^n.n—l 



pdri 
dt n— 1,ti— 1 



G n -2.n-2 



C ^n n—1 j 



dt 2 



t-J-n.n dt n — Ij^ - 1 ' ^ ^n^nj 



k<n-2, 
k = n — 1, 

fc = 71. 



(5.2) 
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Proof 
We start by considering the case k < n — 2 where the first derivative reads 



d 



: / dsi ■ ■ ■ [ 

JO Js k 



\n—k 



k 



I 1 A /"* A ^-^ 

Cr„,fe G„_i,fc + c/ dsi • ■ ■ / ds k — - — — - JLi 

Sfc_i l n fc J- =1 



1 I cosh c(sj — Sj-i) sinhc(i — Sk), (5.3) 



and, from (5.3), the second derivative becomes 



d£^ dt 



(< - **) 



n-fe-1 



+c / dsi • • ■ / ds fc — — — — — — TT cosh c(s., - Sj-i) sinhc(t - s k ) 

Jo Js k _ x (n — « — 1)! 



2-r-G„_i fc — G n -2,k + C G n .k- 

at 



(5.4) 



The expression (5.3), for k = n — 1, is qualitatively different and must be handled carefully. In this case 
we have that 



dt 



(-*n.n — l ^n— l,n- 



pt pt n-1 

/ dsi • • • / ds„_i(t- s n _i) TT coshc(sj - sj_i)sinhc(t- s n _i), (5.5) 



and thus, from (5.5), 
d 2 



, ,o ^"n.n— 1 i.^n— l.n— 1 t C (_r nn — 1 

di^ dt 



+c / dsi • • • / ds„_i 1 | cosh c(sj — Sj-i) sinhc(< — s n _i) 

JO Js„_ 2 j = 1 

-Tj.Gn-l,n-l + ( "T7G n _i :n _i — G n -2,n-2 I +C G n ^ n -\ 



2^r-&n— l.n— 1 — G n _2,n— 2 + C &n,n— 1- 

di 



We omit the derivation of the last formula in (5.2) because it coincides with the result of Lemma 3.1 in 
P. ■ 



The results of Lemma 5.1 permit us to obtain the equation governing the mean- value of the hyperbolic 
distance of the center of mass of the randomly moving splinters produced by the disintegration process. 



Theorem 5.1. The mean-value of (3.2), £ , {coshry cm (i)}, is solution to the non-homogeneus second-order 
linear equation 



-XI 



di 2 



u — c u 



Ac 2 e * 
VA 2 + 2 4 < 



e^ A 



2 4 c 2 1 



Proof 



We begin by successively deriving the expression (5.1) as follows 

, oo n-l 1 , oo .„ , 



dt 



71=1 fc=0 



71=0 



and 



oo n — 1 



1 d 2 



dt 2 



2A^-A 2 . + e--E^E^^ G ^ + e - At E 



71=1 fe=0 



2 k + 1 dt 2 



71=0 



A" d 2 
2™^ 



G r , 



(5.6) 



(5.7) 
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We must now insert all the expressions of Lemma 5.1 into (5.7) by accurately taking into account the 
constrains on n and k. It is convenient to write the last two terms of (5.7) as follows 

oo n— 1 



1 d 



A" d 2 



E A " E ^fc+r ^ G «. fc + E ™ m^ g 



n=l fc=0 

oo n— 2 



l d 



2™di 2 

n=0 
2 - A „ d 2 



E A " E ^FTi ^72 G «> fc + E 



n=2 fc=0 



2" di 2 



+ A d C 



°° A" d 2 d 2 

Ew^ G ™>«+h72 G °>o 



2" di 2 

n-2 
= E A " E ^fc+1 



ra=l 

Co n— 2 



n=2 fc=0 

Z^ 2" 
n=2 

Z_^ On 



2-r-G„_i.fc — G n -2,k + C G n ,k 

di 



d_ 
'di' 



2 — G n _i.„_i — G n _2,n-2 + c Ci„. n _i 



+ A — Gq,0 + X c G^o 



^n— l,n— 1 i C ^n,n 



c 2 G ,o- 



(5.8) 



By regrouping the above terms we notice that (5.8) can take the form 

oo n— 1 1 oo , n 

c 2 ^ a" 22 ^pi G »,fe + E ^ G «^ 



,n=l fe=0 

-A 2 



n=0 



*4 

di 



oo n — 1 ^ oo . n 

E A " E ^fc+i G '^ fc + E ^ G «>" 



oo n— 1 



X A ™ X ^fc+T G «,fc + X ^7 G " 



.n=l fe=0 
d 



n=0 



2 k 
,n=l k=0 

A^A n d 



X 



2 ^ 2™ di 

n=0 



G,, 



n=0 

A 2 °° A" 

22 Z^ 2" 

n=0 



cV'u + 2X^[e xt u] - \ 2 e xt u - J V J- £ 
di L J 2 ^ 2™ di 

n=0 



G n . 



A 2 °° A™ 

22 X 2" G "'"' 

n=0 



(5.9) 



Multiplying by e the right hand side of (5.9), we have that 

oo n— 1 



-A/ 



n=l fe=0 



2 fe +! di 2 ' 



E An E^n^ G «^' 



A/ 



<^2™di 2 n ' n 

n=0 



c 2 u + 2X 2 u + 2\±u - \ 2 u - V A4 V J-f G n „ + ^e- xt V J-G„ „. 

df 9 Z^ On J + ' 92 Z^ On ™' n 



n=0 



n=0 



(5.10) 



By inserting result (5.10) into (5.7) we finally obtain 

A 



—u~c 2 u = --e~ xt , 

at 2 2 <^ 2™ di 

n=0 



^ A" d 



+ * e 



-At 



\ ft 
Z^ ^n" G ™. n ' 



n=0 



(5.11) 



The series e - ^ Sn'Lo W^" n ^ n nas been studied in [T] and represents the mean hyperbolic distance of the 
particle moving in H^ and changing direction at each Poisson event when the rate of the Poisson process 
is A/2. We have that 



A 



Z^ On G ™>" ~ 



e * 



n=0 



1 



A 



\/A 2 + 2 4 c 2 



Since 



-A* ^ A7> d_ 
6 2 Z^ 2 »di "'" 

n=0 



d 
di 



e 2 



,WA 2 +2-ic2 



°° A™ 

Z^ 9^ G "> n 



1 - 



A 



\/A 2 + 2 4 c 2 



,-|VA 2 +2i c 2 



A _ At ■r-^v A™ 

2 / ^ 2™ 

n=0 
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the right-hand side of (5.11 1 can be rewritten as 

A 2 



3G 



2 e 2^2«dt "'" 

A _ xt [ A _ xt 



2 2 



xt \~~^ ™ /~i A _ At d 



n=0 



dt 






A 



VA 2 + 2 4 c 2 _xt 
2^ e 22 



x/A 2 + 2 4 c 2 

A 
\/A 2 + 2 4 c 2 



e £VXS+2*c= 



1- 



A 



VA 2 + 2 4 c 2 



e -^v / A^+2^? 



e iWWW_ (1 _ 



A 



VF+2V 



e -^v^+2^ 



2 4 



_ 3At 

e ~^ r 



A 



VA 2 + 2 4 c 2 



A\/A 2 + 2 4 c 2 _3A* 
-e ~^ r 



2 4 



-e 2^ e 2J 



1 



* VA 2 +2 4 c 2 



A 



A 



£2' 



i v / A 2 +2 4 c 2 



VA 2 + 2 4 c 2 

AC _3At _ _f_ ,, \ J _■)!,. J 

VA 2 + 2 4 c 2 



VA 2 + 2 4 c 2 
1- 



* VA 2 +2 4 c 2 



A 



--r -^-VA 2 +2 4 c 2 



VA 2 + 2 4 c 2 
Ac 2 
VA 2 + 2 4 c 2 ' 



-WA 2 +2 4 c 2 



This lets equation (5.6) appear 



The expression (4.1 ) of the mean- value -E{cosh77 cm (£)} derived in Section 4 satisfies the non-homogeneous 
second-order linear equation (5.6 1 as shown in the next theorem. This proves that the mean hyperbolic 



distance (4.1) is obtained by two different and independent methods. 



Theorem 5.2. The mean hyperbolic distance (4-1) is the solution to the following Cauchy problem, 

j e -^VA 2 +2 4 c 2 _ e ^VAH2Vl 



■prU — C U 



Ac 2 e"^ 



dt 2 " " ™ VA 2 +2 4 c 2 

u(0) = 1, 



d/ 



'(*)U = °- 



Proof 



In order to perform the calculations it is convenient to write (4.1 ) as 



E{cosh Vcm (t)} = KA(t)B(t) + _ A + 2 f , e ct 



A -2c 



2(A + 3c) 2(A-3c) 



where 



K 



r%/A 2 +2 4 C 2 



rVA 2 +2 4 c 2 



x/A 2 + 2 4 c 2 ' 
and B(t) has derivatives 



A{t) = e-^ A \ B(t) = 



5A + 3\/A 2 + 2 4 c 2 5A - 3\/A 2 + 2 4 c 2 ' 



B'(t) 



VA 2 + 2 4 c 2 
2^ 



--WA 2 +2 4 c 2 

e 2 1 



-WA 2 +2 4 c 2 

e^ v 



5A + 3\/A 2 + 2 4 c 2 5A - 3%/A 2 + 2 4 c 2 



A 2 + 2 4 c 2 
, B"{t) = d-±±£-B(ik). 
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Therefore 
-— E{cosh i] cm (t)} - c 2 E {cosh r/ cm (i)} = — — 



r\t 



2 2 Ac 2 e~» At J 5A 
VA 2 + 2 4 c 2 |2 2 

+ d 2 2 

Xc 2 e-^ xt 
VA 2 + 2 4 c 2 



e _^ VA . +2 4 c 2 



r B - 35' 



5A + 3VA 2 + 2 4 c 2 5A - 3VA 2 + 2 4 c 2 



_£VX*+2*c* 



f »^v'A 2 +2 4 c 2 



5A + 3VA 2 + 2 4 c 2 5A - 3VA 2 + 2 4 c 2 



e -|VA 2 +2i C 2 

sA+sv'pn 1 ^ 



fVA 2 +2 4 c 2 

(5A + 3 V / A 2 + 2 4 c 2 ) + — / x9 ^ (-5A + 3VA 2 + 2 4 c 2 ) 



5A - 3VA 2 + 2 4 c 2 



X ° 2e w * [ e -|VA 2 +2 4 c 2 _ p fVA 2 +2 4 c 2 



VA 2 + 2 4 ( 



We can easily check that -E{cosh77 cm (0)} = 1 and we have that ^E{cosh.r] cm (t)}\._ = 



since 



— E{coshr) cm (t)} = KA 



dt 



XL 



„ 3 

2c 2 e & 

VX 2 + 2 4 c 

A + 2c 



-%*+* 



e ^VA 2 +2 4 c 2 

5A- 37^+1^ 



A + 2c 
: 2(A + 3c) 



e ct -c 



A -2c 
2(A-3c)' 



„-ci 



(3A - \/\ 2 + 2 4 c 2 ) - 



e 22 



5A + 3^A 2 + 2 4 c 2 



(3A + \/\ 2 + 2 4 c 2 ) 



A -2c 



'2(A + 3c) 2(A-3c) 



-e ct - c 



Remark 5.1. From Theorem |5.2| an interesting relationship between the hyperbolic mean distance of 
the center of mass and its acceleration can be extracted 



d 2 ^ 2 

' E{cosh r] cm (t)} - — ^E{cosh Tjcmit)} = 



-\t 



t 



dt 2 ' 



VA 2 + 2 4 c 



sinh-r\/A 2 + 2 4 c 2 
2 Z 



(5.12) 



For large values of t, if VA 2 + 2 4 c 2 — 3A > 0, the mean distance and its second derivative tend to coincide 
while in the opposite case they tend to diverge. In view of Remark |4.2| this m eans that there is a different 
rate of growth for the mean distance and its acceleration only if VA 2 + 2 4 c 2 — 3 A < 0. 



6 Mean hyperbolic distance of the k-th splinter 

If N(t) = n> k, the fo-th splinter, at time t, is located at the hyperbolic distance %(£) from O equal to 



fc+i 



coshr7fc(i) = IT COShc(-5j — Sj-i), < k < n, 

3 = 1 



(6.1) 



The expression (6.1) refers to sample paths of particles which change direction until the k-th Poisson 



event and then remain on the same geodesic line until time t. Clearly 770(f) represents the distance of the 
particle which never changed direction and rj n (t) is the distance of the particle that changed direction at 
all Poisson events. 



The mean- value of (6.1 1 becomes 



E{coshrj k (t)I{ N ( t )>k}} = ^2 E{coshT] k (t)I {N(t)=n} } 

n--k 
00 

= J2 E{cosh rj k {t)\N(t) = n}Pr{N{t) = n}. 



(6.2) 



n=k 
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.5 



g=0.05, t=60, N(t)=300 



c-0.05, 1=50. N(t)=300 



c*0.0S, t=50, N(t)r300 



Figure 5: The position of the splinters in the hyperbolic half-plane -Hj at time t = 50, when c = 0.05 and 
N(t) = 300, are drawn. In (a) each splinter chooses the clockwise direction, in (b) each splinter chooses the 
counterclockwise direction and in (c) the clockwise and counterclockwise directions are alternatively chosen. 



We now evaluate the Laplace transform of (6.2 1. By letting 7 = A + \i and by considering (3.4) we have 
that 
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(6.3) 



where f k (s) is the inverse Laplace transform of the function 
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If fc = we have /o(s) = S(s). When k — 1, 2, ... , recalling that L \e QU l+(u) = , _ 1 p+1 with s > a 
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and j = 0, 1, 2 ... , we obtain the inverse Laplace transform 
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By inverting the Laplace transform (6.3) we obtain for k = that E{coshr] (t)Ii N u\ > Q'i } = coshci, while 
for fc = 1, 2, . . . we have 
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where #(s; k, A) = - — °^4 , Y r ^ ~ Beta(r, k — r) and (21^ — 1) € (— 1, 1). If we adopt the convention 

that Yo.fc = 1 an£ i ^fe.fc = —1 it is possible to simplify the last expression rewriting it as 
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where h(k,C,s) — J2r=a { r )^ Y r fc {e cs ' 2y '' fc_1 ^}. The last result suggests regarding the mean hyperbolic 
distance of a particle generated at the fc-th Poisson event as the mean hyperbolic distance of a particle 
moving on the first geodesic line and stopping at a randomly distributed time; the law of the stopping 
time is a suitable combination of a Beta and Gamma distributions. 



Remark 6.1. In some particular cases it is possible to check easily that formula (6.4 1 gives exactly 
the same expression of the mean hyperbolic distance of the k-th particle as that obtained starting from 
formula (6.2). In particular for k = 1 we have 
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For k = 2 we obtain 

E{coshr]2(t)I {N(t) >2}} 

- — / coshc(t-s) e cs +2£y 12 {e cs ( 2yi ' 2 - 1 »} + e- cs 5 (s;2,A)ds 
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- -y cosh c(t-a) cosh cs ds + - J coshc(t - s) E Yl ..{e 81 ^''- 1 '} p d.s 
where Yi 2 ~ Beta(l, 1). 
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